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GAMMA-RAY IMAGING WITH THE CODED MASK IBIS TELESCOPE 
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ABSTRACT 



The IBIS telescope onboard INTEGRAL, the ESA 
gamma-ray space mission to be launched in 2002, is 
a soft gamma-ray (20 keV - 10 MeV) device based 
on a coded aperture imaging system. We describe 
here basic concepts of coded masks, the imaging sys- 
tem of the IBIS telescope, and the standard data 
analysis procedures to reconstruct sky images. This 
analysis includes, for both the low-energy detector 
layer (ISGRI) and the high energy layer (PICSIT), 
iterative procedures which decode recorded shadow- 
grams, search for and locate sources, clean for sec- 
ondary lobes, and then rotate and compose sky im- 
ages. These procedures will be implemented in the 
Quick Look and Standard Analysis of the INTE- 
GRAL Science Data Center (ISDC) as IBIS Instru- 
ment Specific Software. 
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1. CODED MASK IMAGING 



Coded aperture systems nowadays find their major 
application in high energy astronomy, and in partic- 
ular in the hard X-ray (3-30 keV) and soft 7-ray (30 
keV - 20 McV) domains where conventional focusing 
techniques are difficult to implement and where the 
high and variable background limits the performance 
of standard on/off monitoring techniques (see review 
by Carofi et al. 1987). 

In coded aperture telescopes source radiation is spa- 
tially modulated by a mask of opaque and transpar- 
ent elements before being recorded by a position sen- 
sitive detector, allowing simultaneous measurement 
of source and background fiuxes. Reconstruction of 
the sky image is generally based on a correlation pro- 
cedure between the recorded image and a decoding 
array derived from the mask pattern. Essential re- 
quirements for such systems is that the mask pattern 
must allow each source in the field of view (FOV) 
to cast a unique shadowgram on the detector and 



that projected shadowgrams must be, as much as 
possible, different to better differentiate the sources. 
Assuming a perfect detector (infinite spatial resolu- 
tion), the angular resolution of such a system is then 
defined by the angle subtended by one hole at the 
detector. The sensitive area depends on the num- 
ber of all transparent elements of the mask viewed 
by the detector. So, reducing hole size or increasing 
mask-detector distance while increasing accordingly 
the number of holes improves the angular resolution 
without loss of sensitivity. A large field of view can 
also be obtained since the source radiation is modu- 
lated for sources which are within an angle from the 
axis given by the arctg of the mask plus detector di- 
mensions divided by the mask-detector distance. To 
optimize the sensitive area of the detector, masks of 
dimensions greater of or equal to the detector dimen- 
sions are employed. Two kind of FOV are defined. 
The fully coded (FC) FOV for which all source ra- 
diation directed towards the detector plane is mod- 
ulated by the mask and the Partially Coded (PC) 
FOV for which only a fraction of this source radia- 
tion is modulated, while the rest, if detected, cannot 
be distinguished from the background. If holes are 
uniformly distributed the sensitivity is constant in 
the FCFOV and decreasing in the PCFOV. 



2. OPTIMUM CODED APERTURE SYSTEMS 



The System Point Spread Function (SPSF) of coded 
aperture telescopes, i.e. the final imaging response to 
a point source after reconstruction, depends critically 
on the mask pattern. Representing the mask with an 
array M of 1 (open elements) and (opaque ones), 
the detector array D will be given by the convolu- 
tion of the sky image S by M plus an unmodulated 
background array term B, 

Si^M + B 

Suppose to find a special array M for which exists a 
correlation inverse G such that M -k G = (5-function, 
then we can reconstruct the sky by 

S'^Di.G = Si.M*G + Bi.G 
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= Si.6 + BirG = S + Bi.G 

and S' differs from S only by the B -kG term, which 
for a flat array i? is a constant level which can be 
measm'ed and removed. Such special mask patterns, 
including those called uniformly redundant arrays 
(URA), were found in the 70s, and they allow the 
reconstructed image to be free of secondary lobes 
(Fenimore & Cannon 1978). Most of these patterns 
are built using binary sets called cyclic different sets 
which have the remarkable property that their cyclic 
autocorrelation gives a delta function. The decoding 
array G = 2M - 1 (i.e., G=+l for M=l and G=-l 
for M=0) is then a correlation inverse. 
To have a sidelobe-free response a source must be 
able to cast on the detector a whole basic pattern 
(fully coded source) . To make use of all the detector 
area and to allow more than one source to be fully 
coded, the mask basic pattern is normally taken of 
the same size and shape of the detector and the to- 
tal mask made by a cyclic repetition (< 2 x 2 for 
rectangular mask) of the basic pattern. For such op- 
timum systems a FCFOV source will always project 
a cyclically shifted version of the basic pattern and 
correlating the detector image with the G decoding 
array will provide sidelobe-free peak with position- 
invariant shape at the source position. 
A source in the FCFOV will instead cast an incom- 
plete pattern and its contribution cannot be a-priori 
subtracted and will produce secondary lobes (coding 
noise). On the other hand the modulated radiation 
from PC sources can be reconstructed by extending 
with a proper normalization the correlation proce- 
dure to the FCFOV (§4). The complete field of view 
of the telescope (FOV) is therefore composed by the 
central FCFOV of constant sensitivity and optimum 
image properties (position-invariant and flat side- 
lobes SPSF) surrounded by the FCFOV of decreas- 
ing sensitivity and non perfect SFSF. A source out- 
side the FOV simply contributes to the background 
level. 

These masks also minimize the statistical errors of 
the reconstructed peaks. Since V = G^ * D — YjD 
the variance associated with each reconstructed sky 
image pixel is constant in the FCFOV and equal to 
total counts recorded by the detector, therefore the 
source signal to noise is simply 



S/N = Cs/VCs + Cb 

where Cs and Cb are source and background counts. 
These masks also have nearly equal number of trans- 
parent and opaque elements and therefore offer mini- 
mum statistical error in condition of high background 
(typical of the 7-ray domain) . However the sensitiv- 
ity also depends on the detector spatial resolution 
and an imaging efficiency factor must be applied to 
this maximum S/N to account for this effect. 



3. SYSTEM POINT SPREAD FUNCTION 



To perform discrete operations the counts are binned 
in detector pixels and, to avoid loss due to coarse 



sampling, pixels are of much smaller size than the 
mask elements. The correlation can take the form of 
fine cross- correlation, for which the array G is itself 
divided in flner el ements with the same samp ling and 
then correlated ( Fenimore fc Cannon 1981 ). The 
SPSF is then given by a pyramidal function whose 
width (FWHM) is 1 mask element. However the spa- 
tial resolution of real detector is not infinite and this 
induces an intrinsic loss in peak reconstruction and 
makes fine sampling (sampling-pixels small) useless. 
In this case, and in particular when detector resolu- 
tion is not negligible with respect to the mask ele- 
ment size, it can be shown that the best point-source 
signal to noise is obtained by convolving the detec- 
tor image with the G array convolved by the detector 
point spread function. For a pixellated detector the 
blurring function is just a block function of width 1 
pixel. Convolving G with this function and perform- 
ing the correlation we obtain a final SPSF which is 
the convolution of 2 pyramidal functions one of width 
1 mask element and the other of width 1 pixel. 
Due to finite pixel size, the peak of the SPSF will be 
reconstructed in the best case with an average effi - 
ciency given by the relation (1—3^) ( Skinner 1995 ), 
where R is the ratio of the mask element size to pixel 
size. This gives the average loss in sensitivity when 
the flux is estimated by fitting the image with the 
SPSF function at fixed source position. On the other 
hands when performing the SPSF deconvolution with 
given sampling the reconstructed peak height will be 
even lower, due to the fact that the source will not 
be at the center of the sampled pixel. In average 
this (imaging) efficiency will be given by (1 — jj^)^- 
In the case of continuous detectors, like the SIGMA 
telescope, the detector blurring can be described by 
a bi-dimcnsional Gaussian, and the SPSF by an an- 
alytical function (e.g. appendix A in Bouchet et al. 
2001), which depends on the detector spatial resolu- 
tion (width of the Gaussian). To precisely evaluate 
source parameters and their errors we compare re- 
constructed image sectors with the system PSF by 
means of the chi-square fitting technique. 



4. DECONVOLUTION AND ANALYSIS 



Discrete cross-correlation to compute sky and vari- 
ance images can be written 



Eg 



i+k,j+iDki ; 



kl 



kl 



where Poisson statistics was assumed. This standard 
deconvolution in FCFOV can be extended in the FC- 
FOV by extending the correlation of the decoding ar- 
ray G with the detector array D in a non-cyclic form, 
padding G with elements. Since only the detector 
section modulated by the PC source is used to recon- 
struct the signal, the statistical error at the source 
position and significance of the ghost peaks are mini- 
mized. However to ensure a fiat image in the absence 
of sources, detector pixels which for a given sky po- 
sition correspond to mask opaque elements must be 
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balanced, before subtraction, with the factor b — ^ 

. . . 

where is the number of pixels corresponding to 

transparent elements and n_ to opaque ones for that 

given sky position. This can be written 



where the decoding arrays are obtained from the 
mask M by G+ = M and = 1 - M, then padded 
with O's outside mask region, and where the sum 
is performed over all detector elements. In the FC- 
FOV we obtain the same result of the standard cross- 
correlation. To consider effects such as satellite drift 
corrections (see Goldwurm 1995), dead areas or other 
specific conditions, a weighting array W is used to 
weigh properly the detector array before correlating 
it with the G arrays. The balance array is 



Y.kGi 



EkG:, 



+k,j+l 



The variance, which is not constant outside the FC- 
FOV, is computed accordingly 



Wki 



B 



2. V 



Dki (' 



G 



i+k,j+l 



Wkl 



since the cross-terms G^G~ vanish. Note however 
that when the weights Wki refer to the same pixel 
in D, the terms G W must be summed before squar- 
ing (see Goldwurm 1995). The varying effective area 
can be calculated by similar formula and used to 
renormalize, after background subtraction to FC- 
FOV count rates. All this can be performed for sam- 
pling finer than 1 pixel per mask element and using 
a G array convolved with detector PSF in order to 
optimize S/N for point sources, with corresponding 
normalizations. This procedure can be carried out 
with a fast algorithm by reducing previous formulae 
to a set of correlations computed by FFT. 
The on-axis SPSF on the whole FOV for an 
optimum system (IBIS/ISGRI configuration with 
MURA mask and R = 2.43, § 5) and PSF deconvolu- 
tion is shown in Fig. ^ Note the peak and flat level 
in the central FCFOV, the secondary lobes (coding 
noise) in the PCFOV and the 8 main ghosts of the 
source peak in the PCFOV located at distances, from 
the source, which are multiple of the basic pattern. 



5. IBIS IMAGING SYSTEM 



The IBIS coded mask imaging system is composed 
by a rephcated Modified URA (MURA) mask of 
tungsten elements (Fig. ^ and 2 pixellated detec- 
tor planes of the same size, ISGRI for the low 
energy band (20-1000 keV) and PICSIT for the 
higher band (1 50 keV - 10 MeV) disposed about 
10 cm below (lUbertini et al 1998[) . The MURA 
( Gottesman fc Fenimore 1989|) are nearly-optimum 




Figure 1. SPSF for the IBIS/ISGRI telescope. 




Figure 2. The 53 x 53 MURA basic pattern. 



masks and previous discussion for such system is 
valid for the IBIS telescope. The detector planes can 
be divided in a regular grid of square pixels where 
each detector element occupies a pixel of the array. 
Essential imaging characteristics and performances 
are reported in Table 0, including value of cfhciency 
for fitting procedure and imaging. 



6. IBIS POINT SOURCE LOCATION ERROR 



The average Point Source Location Error (PSLE) for 
an optimum coded aperture system with a defined 
System Point Spread Function depends on source sig- 
nal to noise ratio (S/N) as following 

PSLE - ^ 



R (S/N) 



This can be demonstrated by computing, through 
the covariance matrix, the errors on the parame- 
ters (source position and intensity and background 
level) derived by the chi-square fitting of a decon- 
volved image sector with the SPSF. We calculated 
this error, averaged over a large number of uniformly 
distributed source positions, for an optimum coded 
aperture system where mask elements are an inte- 
ger number of detector pixels for different values of 
mask/pixel sizes and S/Ns. In Fig. ^ are reported 
these theoretical location errors for mask to pixel ra- 
tios R = 2 and 3, which are respectively upper and 
lower limits for the IBIS/ISGRI configuration which 
has an intermediate ratio of i? = (11.2/4.6) = 2.43. 
In the same plot are reported location errors ob- 
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BIS/ISGRI POSITION ERROR vs. S/N 



Table 1. IBIS Imaging System and Performances 



MURA basic; total mask 
Mask element size 
ISGRI-Mask distance 
ISGRI active pixels 
ISGRI pix size 
ISGRI pixel pitch 
ISGRI pixels 
PICSIT-Mask distance 
PICSIT Active Pixels 
PICSIT pix size 
PICSIT pixel pitch 
PICSIT total pixels 

ISGRI effic. (Fit-Imag) 
ISGRI ang.res.(FWHM) 
ISGRI pixel angle 
PICSIT effic. (Fit-Imag) 
PICSIT an.res.(FWHM) 
PICSIT pixel angle 
EXFOV (0 % sens.) 
FCFOV (100 % sens.) 



53 X 53 ; 95 x 95 

11.2 X 11.2 X 16 mm^ 

3200 mm (top — top) 

128 X 128 

4. X 4. X 2. mm^ 

4.6 X 4.6 mm^ 

130 X 134 

3300 mm {top — top) 

64 X 64 

8.66 X 8.66 X 30 mm^ 



9.2 X 

65 X 

0.86 - 
12.0' 
5.0' 
0.73 - 
11.7' 
9.6' 
30.6° 
8.3 



9.2 

67 

0.81 



0.63 



mm 



X 31.0° 
X 8.6° 



tained from analysis of simulated observations with 
IBIS/ISGRI system of an on-axis point source. 200 
ISGRI images for a point-like nearly-on-axis source 
of given S/N (the S/N is here reduced by imaging ef- 
ficiency as defined in §3) are simulated, deconvolved 
and analyzed to search for the most significant peak 
which is then fitted to the analytical SPSF to derive 
source position. Background (uniform distribution) 
was kept fixed to a value of 1000 cts and S/N was 
made varying by reducing input source counts from 
w 2000 to « 500. The standard deviation of the 
observed offsets between the best fit positions and 
the input positions represent then the 1 a error in 1 
parameter of the estimated position. These results 
show that the estimated location error, as expected, 
varies with S /N and R as predicted and can therefore 
be for bright sources a small fraction of the angular 
resolution, even for small values of the ratio mask el- 
ement size to pixel size. The IBIS/ISGRI telescope, 
assuming no error in pointing axis reconstruction or 
other systematic effects, can locate a 30 cr point-like 
source at better than ± 30". Absolute error in at- 
titude reconstruction for INTEGRAL is expected to 
be < 20". 



7. ITERATIVE IMAGE RECONSTRUCTION 



In a standard analysis, IBIS events or histograms 
are binned in detector images, which are then cor- 
rected for detecto r and background non-uniformity 
(Goldwurm 1995) and then processed by an itera- 
tivc algorithm which decodes, cleans and composes 
sky images. For each detector image a sky image 
and its variance are obtained using the deconvolu- 
tion procedure of §4, and then iteratively searched 
for sources and cleaned of the source side lobes. In 



SlmuL {Bkg= 1 000 cts) 
Predic. Upper Lim. 
Predic. Lower Lim. 




20 40 
Signal to Noise Ratio 



Figure 3. Position errors ( at QCfVo confidence level 
in 2 parameters) for the IBIS/ISGRI telescope from 
simulations of a nearly on-axis source at different 
S/N compared to predicted theoretical limits. Solid 
line is a function (S/N)^^ normalized at S/N=9.5. 



this iterative process the source peaks are fitted with 
the proper system PSF and finely located. Then the 
source contribution to the image is modeled in detail 
and subtracted. The images are rotated, projected 
and summed after beeing weighted with the variance, 
and then searched for further contributions. Note 
that the search for significant excesses must be per- 
formed taking into account that these are correlation 
images and the critical level at which an unknown ex- 
cess is significant incre ases typically from 3 cr to 5 - 
6 a ( [Caroli et al. 198*7 ). This procedure was used in 
ulat 



simulations to assess capability of IBIS /ISGRI tele- 
scope for a number of specific scientific cases (see e.g. 
Goldwurm et al. 1999, 2000, and Goldoni et al. 2000 
these proceedings) , and it is now being implemented 
in software modules to be integrated into the ISDC 
scientific analysis pipelines. 
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